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INTRODUCTION 


Based  on  a  detailed  study  reported  in  Reference  1,  the  best 
material  model  for  gun  tubes  under  high  pressure  operation  is  an 
elastic-plastic  material  which  obeys  the  Mises'  yield  criterion  and 
the  Prandtl-Reuss  incremental  stress-strain  relation.  A  literature 
survey  indicates  that  no  closed  form  solution  exists  even  for  the 
axisymmetric  plane-strain  problems.  And,  in  such  situations,  one  has 
to  rely  on  numerical  methods.  Both  the  finite-difference  method^’ ^ 
and  the  finite-element  method4  have  been  used  to  solve  the  elasto- 
plastic  problems  considered  here.  The  finite-element  method  is  more 
powerful  and  can  be  used  to  solve  more  general  elastoplastic  problems. 
Since  the  displacement  function  is  assumed  and  the  programming  is 
complicated,  the  accuracy  of  the  finite-element  approach  has  to  be 
verified.  This  is  usually  done  by  comparing  more  rigorous  solutions 
to  simpler  problems.  For  the  problem  considered  here,  rigorous 

solutions  based  on  the  finite-difference  method  were  obtained  by 

2  ? 

Hodge  and  White  for  ideally-plastic  materials  and  by  Chu  for  strain¬ 
hardening  materials. 


^Davidson,  T.  E.,  and  Kendall,  D.  P.,  "The  Design  of  Pressure  Vessels 
For  Very  High  Pressure  Operation,"  Watervliet  Arsenal  Report  WVT-6917. 
Also  in  Mechanical  Behavior  of  Materials  Under  Pressure  (edited  by 
Pugh,  H.  L.  D.),  Elsevier  Co.,  1970,  Chapter  2. 

^ White,  P.  G.  and  White,  G.  N.,  "A  Quantitative  Comparison  of  Flow  and 
Deformation  Theories  of  Plasticity,"  J.  Appl.  Mech.,  Vol.  17,  1950, 
pp.  180-184. 

Chu,  S.  C.,  "A  More  Rational  Approach  to  the  Problem  of  an  Elasto- 
Plastic  Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol.  294, 
1972,  pp.  57-65. 

4Chen,  P.  C.  T„,  "The  Finite  Element  Analysis  of  Elastic-Plastic  Thick- 
Walled  Tubes,"  Proceedings  of  Army  Symposium  on  Solid  Mechanics,  1972, 
The  Role  of  Mechanics  in  Design-Ballistic  Problems,  pp.  243-253. 
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In  the  present  paper,  a  new  finite-difference  approach  is  developed 
for  solving  the  axisymmetric  plane  strain  problems  subjected  to  internal 
or  external  pressure  beyond  the  elastic  limit.  An  incremental  approach 
is  used,  and  the  numerical  scheme  is  stable  for  ideally-plastic  as  well 
as  strain-hardening  materials.  The  approach  is  simpler  than  others, 
yet  very  accurate  solutions  can  be  obtained  by  reducing  the  grid  sizes 
and  load  increments. 

BASIC  EQUATIONS 

Assuming  small  strain  and  no  body  forces  in  the  axisymmetric  state 
of  plane  strain,  the  radial  and  tangential  stresses,  ar  and  a0,  must 
satisfy  the  equilibrium  equation, 

r(3ar/3r)  =  aQ  -  ar  ;  (1) 

and  the  corresponding  strains,  er  and  e0,  are  given  in  terms  of  the 
radial  displacement,  u,  by 

er  =  3u/3r  ,  eQ  =  u/r  .  (2) 

It  follows  that  the  strains  must  satisfy  the  equation  of  compatibility 

r (3e0/3r)  <=  er  -  eQ  .  (3) 

Whereas  the  differential  equations  (1),  (2)  and  (3)  hold  throughout  the 
tube  regardless  of  the  material  properties,  the  constitution  equations 
assume  various  forms  according  to  the  adopted  form  of  yield  function, 
hardening  rule,  total  or  incremental  theory  of  plasticity.  In  the 
present  paper,  the  material  is  assumed  to  be  elastic-plastic,  obeying 
the  Mises'  yield  criterion,  the  Prandtl-Reuss  flow  theory  and  the 
isotropic  hardening  law.  The  complete  stress-strain  relations  are:5 

“i 11,  R.,  Mathematical  Theory  of  Plasticity,  Oxford  University  Press. 
1950. 


2 


(4) 


det '  =  dai,/2G  +  (3/2)ai'da/(aH') 

da  _>  0  for  i  =  r,0,  z 

de  =  E“1(l-2v)da 
m  v  J  m 

where  E,  v  Young’s  modulus,  Poisson's  respectively, 

2G  =  E/(l+v) 

em  =  (£r+ee+ezV3  ,  , 


a  =  (a  +aQ+a  )/3  ,  a. '  =  a. -a  , 

m  v  r  6  z'-  *  i  x  m  * 


a  =  (1//2)  [(ar-aQ)2  +  (a0-az)2  +  (az-ap2]1/2  >;  aQ  , 
and  aQ  is  the  yield  stress  in  simple  tension  or  compression.  For  a 
strain  hardening  material,  H*  is  the  slope  of  the  effective  stress/ 
plastic  strain  curve 


(5) 


(6) 


a  =  H(JdeP)  .  (7) 

For  an  ideally-plastic  material  (H'=0),  the  quantity  (3/2)da/(aH')  is 
be  replaced  by  dX,  a  positive  factor  of  proportionality.  When  O  <  oQ 
or  da  <  0,  the  state  of  stress  is  elastic  and  the  second  term  in  equa¬ 
tion  (4)  disappears.  Following  Yamada  et  al,^  equations  (4)  and  (5) 
can  be  rewritten  in  an  incremental  form 

dai  =  dijdej  for  =  r,0»z 


and 


dij/2G  =  v/(l-2v)  +  -  a.'a.'/S 


(8) 


^Yamada,  Y.,  Yoshimura,  N.,  and  Sakurni,  T.,  "Plastic  Stress-Strain 
Matrix  and  Its  Application  for  the  Solution  of  Elastic-Plastic  Prob¬ 
lems  by  the  Finite  Element  Method,"  Int.  J.  Mech.  Sci.,  Vol.  10,  1968, 
pp.  343-354. 
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where 


S  =  |  [1  +  jH'/G)a2  , 


(9) 


and  6^^  is  the  Kronecker  delta. 


This  form  was  used  in  the  finite-element  formulation  for  solving 

4 

elastic-plastic  thick-walled  tube  problems.  In  the  following  section, 
the  incremental  stress-strain  matrix  will  be  used  in  the  finite  differ¬ 
ence  formulation. 

FINITE-DIFFERENCE  FORMULATION 

Consider  a  thick-walled  cylinder  of  inner  radius  a  and  external 
radius  b.  The  tube  is  subjected  to  inner  pressure  p  and/or  external 
pressure  q.  The  elastic  solution  for  this  problem  is  well-known  and 
the  pressure  p*  or  q*  required  to  cause  initial  yielding  can  be  deter¬ 
mined  by  using  the  Mises'  yield  criterion.  For  pressure  beyond  the 
elastic  limit,  an  incremental  approach  of  the  finite-difference  formu¬ 
lation  is  used.  The  analysis  starts  with  the  applied  pressure  p  or  q 
and  the  loading  path  is  divided  into  m  increments  with 


Ap  =  (p-p*)/m  ,  Aq  =  (q-q*)/m  . 


(10) 


The  cross  section  of  the  tube  is  divided  into  n  rings  with 


(11) 


^Chen,  P.  C.  T.,  "The  Finite  Element  Analysis  of  Elastic-Plastic  Thick- 
Walled  Tubes,"  Proceedings  of  Army  Symposium  on  Solid  Mechanics,  1972, 
The  Role  of  Mechanics  in  Design-Ballistic  Problems,  pp.  243-253. 
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where  p  is  the  radius  of  the  elastic-plastic  interface.  At  the  begin¬ 
ning  of  each  increment  of  loading,  the  distribution  of  displacements, 
strains  and  stresses  are  assumed  to  be  known  and  we  want  to  determine 
Au,  Ae^.,  Ae0,  Aar,  Aa0,  Aaz  at  all  grid  points.  Since  the  incremental 
stresses  are  related  to  the  incremental  strains  by  the  incremental 
form  (Eq.  (8))  and  Au  =  rAe0,  there  exists  only  two  unknowns  at  each 
station  that  have  to  be  determined  for  each  increment  of  loading.  The 
unknown  variables  in  the  present  formulation  are  (Ae0)i,  (Ae^^  for 
i  =  l,2,...n,n+l. 

The  equation  of  equilibrium  (1)  and  the  equation  of  compatibility 
(3)  are  valid  for  both  the  elastic  and  the  plastic  regions  of  a  thick- 
walled  tube.  The  finite-difference  forms  of  these  two  equations  at 
i  =  l,...,n  are  given  in  Reference  3  by 

trur2ri)<tor)i  -  fri.rri)(AVi  *  ri<iar)i.i 

‘  [Vfri)(VVi  -  V'Vi.l  "  ‘Vi1  tl2) 

for  the  equation  of  equilibrium,  and 

cri.r2ri>«Vi  -  (ri.rV tAVi  *  ri<AVi*i 

"  (ri*l"ri)(Ve6)i  -  VtePi.i  -  Ce6)i]  (13) 

for  the  equation  of  compatibility. 
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Plastic  Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol. 
294,  1972,  pp.  57-65. 
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With  the  aid  of  the  incremental  stress-strain  relations  (Eq.  (8)), 
equation  (12)  can  be  rewritten  as 

^ri+r2ri^d12^i  +  ^~ri+l+ri^d22^i^Ae0^i 

+  [(ri+r2ri)(dlPi  +  ^i+l^^PJ^Pi 

+  ri^d12^  i+1  ^Ae0^  i+1  +  ^^ll^i+l^Pi+l 
=  Cri+rri)  (a0-ar)i  -  r±  [  Cor) i+1  -  (a^]  .  (14) 

The  boundary  conditions  for  the  problem  are 

Aar(a,t)  =  -Ap  ,  Aar(b,t)  =  -Aq  .  (15) 

Using  the  incremental  relations  (Eq.  (8)),  we  rewrite  equation  (15)  as 

(d12^i^^0-^i  +  (dll)  1  (Aer)  i  =  ~^p  >  (16) 

and 

(^n+l^n+l  +  ^dll^n+l (Aer^ n+1  =  *  (17) 

Now  we  can  form  a  system  of  2 (n+1)  equations  for  solving  2 (n+1) 

unknowns,  (Aeq)^,  (Aer)^,  for  i  =  1,2, . . . ,n,n+l.  Equations  (16)  and 
(17)  are  taken  as  the  first  and  last  equations,  respectively,  and  the 
other  2n  equations  are  set  up  at  i  =  l,2,...,n  using  (13)  and  (14). 

The  final  system  is  an  unsymmetric  band  matrix  with  the  nonzero  terms 
clustered  about  the  main  diagonal,  two  below  and  one  above.  The 
storage  required  is  8 (n+1)  words  for  this  system  and  at  least  35 (n+1) 

3 

words  for  the  approach.  In  the  computer  program  which  was  developed, 
the  Gaussian  elimination  method  was  used  to  solve  these  equations. 

All  calculations  were  carried  out  on  IBM  360/Model  44  with  double 
precision  to  reduce  round-off  errors. 

^Chu,  S.  C.,  "A  More  Rational  Approach  to  the  Problem  of  an  Elasto-Plas 
tic  Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol.  294, 
1972,  pp.  57-65. 
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NUMERICAL  RESULTS 


The  plane-strain  problems  of  thick-walled  cylinders  subjected  to 
internal  pressure  p  or  external  pressure  q  beyond  the  elastic  limit 
were  solved.  The  elastic-perfectly-plastic  as  well  as  strain-harden¬ 
ing  materials  were  considered.  The  numerical  results  were  based  on 
the  following  parameters:  b/a  =  2,  v  =  0.3,  H'  =  0  or  E/19.  Various 
values  of  m  and  n  were  used  to  test  the  convergence  of  the  numerical 
results.  The  incremental  loadings  were  applied  until  the  fully 
plastic  state  was  reached.  The  values  for  p  or  q  corresponding  to 
this  final  state  were  denoted  by  p**  or  q**.  It  was  found  that  the 
results  for  these  values  converge  by  increasing  m  and/or  n.  For  an 
ideally-plastic  tube  under  internal  pressure  only,  we  have  (n,m,p**/a0) 
=  (20,  1277,  0.81341),  (50,  1000,  0.80455),  (100,  998,  0.80155),  (200, 
920,  0.80029).  Since  the  storage  needed  for  this  approach  is  only  22% 
of  that  in  Reference  3,  much  larger  n  can  be  used  to  yield  better 
results.  Additional  results  shown  in  Figures  1  to  5  were  based  on 
n  =  100  and  m  =  900-1000.  Figure  1  shows  the  relation  between  inter¬ 
nal  pressure  p,  external  pressure  q  and  elastic-plastic  boundary  p  in 
an  elastic-perfectly-plastic  tube.  Figure  2  shows  the  bore  radial 
and  tangential  strains  as  functions  of  internal  pressure  p  in  an 

J  1 

Chu,  S.  C.,  "A  More  Rational  Approach  to  the  Problem  of  an  Elasto- 
Plastic  Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol. 

294,  1972,  pp.  57-65. 
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ideally-plastic  as  well  as  a  strain-hardening  tube.  Figure  3  shows 
the  distributions  of  radial,  tangential  and  axial  stresses  in  an 
ideally-plastic  tube  subjected  to  internal  pressure.  The  dotted 
curves  correspond  to  initial  yielding.  The  solid  and  broken  curves 
correspond  to  the  cases  when  half  and  all  of  the  tube  is  plastic, 
respectively.  For  an  ideally-plastic  tube  subjected  to  external 
pressure  q,  the  results  were  presented  graphically  in  Figures  1,  4, 
and  5. 
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Figure  1.  Internal  pressure  p  or  external  pressure  q  as  functions  of 
elastic-plastic  boundary  p  in  an  ideally-plastic  tube 
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Figure  2.  Bore  radial  and  tangential  strains  as  functions  of  internal  pressure 
(a)  ideally-plastic  tube,  (b)  strain-hardening  tube,  H'  =  e/19 
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Figure  4.  Bore  radial  and  tangential  strains  as  functions  of  external  pressure 
(a)  ideally-plastic  tube,  (b)  strain-hardening  tube,  H'  =  E/19 
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Figure  5.  Distributions  of  radial,  tangential  and  axial  stress  components  in  an  ideally-plastic 
tube  subjected  to  external  pressure  (p/a=  1.0 . .  1.5  _ ,  2.0  —  —  ) 
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